Anisotropic transport in the two-dimensional electron gas in the presence of 

spin-orbit coupling 
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In a two-dimensional electron gas as realized by a semiconductor quantum well, the presence of 
spin-orbit coupling of both the Rashba and Dresselhaus type leads to anisotropic dispersion relations 
and Fermi contours. We study the effect of this anisotropy on the electrical conductivity in the 
presence of fixed impurity scatterers. The conductivity also shows in general an anisotropy which 
can be tuned by varying the Rashba coefficient. This effect provides a method of detecting and 
investigating spin-orbit coupling by measuring spin-unpolarized electrical currents in the diffusive 
regime. Our approach is based on an exact solution of the two-dimensional Boltzmann equation and 
provides also a natural framework for investigating other transport effects including the anomalous 
Hall effect. 



I. INTRODUCTION 

In the recent years the emerging field of spintronics^^^ 
has generated an intense interest in effects of spin- 
orbit interaction in low-dimensional semiconductor het- 
erostructures. For conduction band electrons in zinc- 
blende semiconductors the dominant effects of spin-orbit 
interaction in low-dimensional geometry can be described 
in terms of two effective contributions to the Hamilto- 
nian. On the one hand there is the Rashba spin-orbit 
term"^ which is due to the inversion-asymmetry of the 
confining potential and has the form 



(1) 



where p is the momentum of the electron confined in 
a two-dimensional geometry, and a the vector of Pauli 
matrices. The coefficient a is tunable in strength by 
the external gate perpendicular to the plane of the two- 
dimensional electron gas. The other contribution is the 
Dresselhaus spin-orbit term which is present in semicon- 
ductors lacking bulk inversion symmetry When re- 
stricted to a two-dimensional semiconductor nanostru- 
ture grown along the [001] direction this coupling is of 
the form^'^ 



(2) 



where the coefficient /3 is determined by the semiconduc- 
tor material and the geometry of the sample. 

The interplay of these two types of spin-orbit cou- 
pling has been investigated theoretically with respect 
to several physical phenomena including the spin- 
splitting in zero magnetic field^~^°, spin precession 
and relaxation^"'^'^'^^^^^, contributions to magneto- 
oscillations^*, quantum interference corrections to the 
conductivity-'^^'^'', mesoscopic transport through quan- 
tum dots^^ , and issues of gate operations between quan- 
tum dot spin-qubits^^ . In the present work we point 
out another effect occurring in the presence of both 



Rashba and Dresselhaus spin-orbit coupling, namely an 
anisotropy in the electrical conductivity for diffusive spin- 
unpolarized charge transport in the presence of fixed im- 
purity scatterers. 

One of the most important concepts in the discussion 
of spin-orbit coupling in semiconductors, and in the field 
of spintronics in general, is the spin field-effect transistor 
proposed by Datta and Das^^. This proposal uses the 
Rashba spin-orbit coupling to perform controlled rota- 
tions of spins of electrons passing through an FET-type 
device operating in the ballistic transport regime. Bal- 
listic transport is necessary to avoid randomization of 
the spin state by even spin-independent scatterers which 
would change the electron momentum and therefore also 
the effective field provided by the Rashba term in an un- 
controlled way. The requirement of ballistic transport 
has been so far one of the major obstacles toward the 
practical realization of a spin FET. Recently an alterna- 
tive scenario for a spin FET was proposed which can also 
operate in the nonballistic regime^^. This proposal ex- 
ploits the fact that if the Rashba coefficient is tuned via 
external gates to be equal to the Dresselhaus coefficient, 
a = /?, a new conserved quantity arises which prohibits 
the randomization of the spin. As we shall see below, 
this particular point a — f3 in parameter space will also 
be of special interest in our present study. 

Our description of diffusive two-dimensional transport 
in the presence of spin-orbit interaction and fixed impuri- 
ties is based on an exact solution of the two-dimensional 
Boltzmann equation and provides also a natural frame- 
work for investigating other transport effects including 
the anomalous Hall effect. In particular, our scheme of 
generating exact solutions to the transport equation can 
deal with arbitrary dispersion relations and is not re- 
stricted to isotropic cases. 

Another study on the possible influence of the Dres- 
selhaus spin orbit coupling on the operation of the spin 
FET was performed very recently by Lusakowski, Wro- 
bel, and Dietl^^. These investigations are restricted to 
the ballistic regime, but take into account, in addition 
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to the Hamiltonian (2), also contributions to the Drcs- 
selhaus term being trilinear in the momentum. In the 
present study we wiU mostly neglect these trilinear con- 
tributions, but discuss their possible influence briefly in 
section IV C. In another very recent preprint R. Winkler 
also studied spin transport in the presence of spin-orbit 
interaction steming from both structure inversion asym- 
metry and bulk inversion asymmotry^^. In particular, in 
Refs.^®'^^ the possibilities arising from different growth 
directions for the two-dimensional electron system are ex- 
plored. Yet another recent work dealing with transport 
in the presence of spin-orbit coupling was performed by 
Mishchenko and Halperin who derived the equations of 
motion for the single electron density matrix in Wigner 
representation in a two-dimensional free electron gas^''. 
The authors applied their results to the dynamic con- 
ductivity of the system taking into account, however, the 
Rashba term only such that no anisotropy in conductiv- 
ity occurred. Finally we mention a very recent work by 
Ganichcv et al. who present an experimental method to 
distinguish the effects of Rashba and Dresselhaus spin- 
orbit coupling using the spin galvanic effect^^. 

This paper is organized as follows. In section II we re- 
view the dispersion relations and eigenstates of free elec- 
trons confined in two dimensions in the presence of spin 
orbit coupling of both the Rashba and the Dresselhaus 
type, and we present results for the Fermi contours. In 
section III we present a scheme to generate exact solu- 
tions to the two-dimensional Boltzmann equation that 
underlies our present study. This approximation-free so- 
lution to the scmiclassical transport equation is then ap- 
plied in section IV to the case of free electrons being 
subject to spin-orbit interaction of the above type. We 
close with a summary and discussion of the results in 
section V. 



II. DISPERSION RELATIONS, EIGENSTATES, 
AND FERMI CONTOURS 

We consider the single-particle Hamiltonian for a two- 
dimensional electron system 



Ti- = h Ti-R + Hd 

2m 



(3) 



where m is an effective band mass. The eigenenergies are 
given by 



(fc) = ^ ± \/ {aky + 0k,f + {ak, + f3kyf (4) 



2m 

with eigenstates 

{r\k,±) 



V2 V 

where A is the area of the system and 



(5) 



X{k) = axg{-aky - pk^ + i{akx + I3ky)) . (6) 
The semiclassical particle velocities are given by 
de± (fc) %k 



v±{k) 



ndk 



m 



± 



[a"^ +I3'^)k + 2al3 (cr^fc) 



(7) 



where is a usual Pauli matrix acting on the vector 
k. As a consistency check let us consider the quantum 
mechanical velocity operator 



[HA ■ 



(8) 



Using the above expressions for the eigenstates it is 
straightforward to show that its matrix elements are 
given by 



{k,±\i^k',±)=Srpv±{k) 



(9) 



i.e. the semiclassical velocities v±{k) are, as usual, the 
diagonal elements of the velocity operator. 

Parametrizing wave vectors as fc = k{cosip,sm(p) 
one obtains for positive Fermi energy £/ the following 
parametrization of the Fermi contours: 



ki{^;sf) = T\ {^) {a^+l3^ + 2a/3sin {2<fi)) 



(a2+/32 + 2a/3sin(2<p)). (10) 



Here the double sign corresponds to the above two dis- 
persion branches, and the Fermi wave vector is given by 



ki {ip;ef) = ki (ifief) (cos sin <^) . 



(11) 



At negative Fermi energies the Fermi contours can be- 
come somewhat more complicated. This case corre- 
sponds to rather low electron densities and shall not be 
considered here further. In the following the Fermi en- 
ergy is always assumed to be positive. Prom Eq. (10) one 
finds the electron density n as 



n ■ 



(12) 



If a = or /3 = the dispersions are isotropic 
and Fermi contours are concentric circles. For a ^ 
7^ /? the Fermi contours arc anisotropic which, as we 
shall see below, leads in general to anisotropic trans- 
port properties. Note that the dispersion relations and 
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Fermi contours arc symmetric around the points (p G 
{7r/4, 37r/4, 57r/4, 77r/4}, i.e. these quantities are invari- 
ant under reflections along the (1,1) and (1,-1) direc- 
tion. These directions define the symmetry axes of the 
problem. In particular, for these direction the wave vec- 
tors and particle velocities are collinear. 

The above findings for the Fermi contours are illus- 
trated in Fig. 1 where we show data for typical values 
for the Fermi energy, Drcsselhaus coefficient and effec- 
tive band mass^^~^® at various values for a^^~^'^. If both 
a and /3 are nonzero the Fermi contours are anisotropic 
having the (1, 1) and (1, —1) direction as symmetry axes. 
The case a = ±/3 is particular^^. Here a new conserved 
quantity given by S := [a^' =F o"'')/\/2 arises'^"'', and the 
spin state of the electrons becomes independent of the 
wave vector. For this situation the dispersion relations 
are more conveniently written as (choosing a = +/?) 



where 



^^±^(M) 



(13) 



(14) 



is the distance vector between the centers of the circles to 
points on their circumference. The double sign labeling 
branches in Eq. (13) does not correspond to the one in 
Eq. (4) and is therefore put in parentheses. As seen in 
Fig. 1, right at a = (3 different parts of the dispersion 
branches for a ^ (3 merge to different circles inducing a 
relabeling of branches. 



III. BOLTZMANN THEORY OF ANISOTROPIC 
TRANSPORT IN TWO DIMENSIONS 

The Boltzmann eqiiation for transport in the two- 
dimensional electron gas in the presence of fixed random 
impurities reads^^ 



dt dr^'' dk~\dt 



(15) 



coll 



Here is a band index; in the context of the previous 
section it corresponds to the double sign labeling the two 
dispersion branches. fn{r,k,t) is the usual semiclassical 
distribution function, and the collision term is given by 

dt Aon 

-w{k',fi';k,ii)f^{k) (l - /^,(fc')) ] , (16) 

where w(fc, /i; fc', //) is a transition probability deter- 
mined by the fixed impurities. The semiclassical equa- 
tions of motion read 



: Vf,{k) 



hdk 



%k 



-eE, 



(17) 



where £^(A;) is the dispersion of the band /U, (— e) = — |e| 
is the electron charge, and E is an external electric field in 
the plane of the two-dimensional gas. Assuming a homo- 
geneous system in a stationary state, fn{r,k,t) = /^(fc), 
and elastic scattering fulfilling the microreversibility con- 
dition w{k, fi;k' , iJ,') = w{k' , iJ,';k, IJ.), the Boltzmann 
equation^^ becomes in lowest order in \E\ 

-eE-v^{k)' 



de 

with the scattering operator 

S [U{k) = 
(fk' 



S 



(18) 



E 



(2^)2 



w{k,t,;k',l,') (/^(fc) - U'{k'))] . (19) 



Here is the equilibrium Fermi distribution depending 
only on the energy e, and the derivative in Eq. (18) has 
to be evaluated at e = £^{k). 

Now let 'd{d) be the angle a given vector a forms with 
the direction of E, fulfilling the relations 



E ■ a = Eacos{'&{d)) , 
X E] ■a = Easin{{){a)) . 



(20) 
(21) 



where is the direction perpendicular to the two- 
dimensional (xy)-plane. The form of the transport equa- 
tion (18) suggests to study the action of the scattering op- 
erator (19) on the functions f^{k) = \v^{k)\ cos[i?(i'^(fc))] 

and f^i{k) = \v^{k)\sm['&{v^{k))]. In fact, by insert- 
ing these functions into (19) and expressing the angle 

t?(v(fc')) intermsof ^?(i;^(A?)) and (^?(t;^(fc))-t?(v (^'))) 
via elementary trigonometric relations, it is easy to show 
that 

\v^{k)\cos[^{v^{k))]\ 
\v^ik)\sin[d{Mm J. 

~^)\ f \Mk)\cos[^iMk))] 
\v^ik)\sin[{){v^ik))] 



(22) 



with 



rjlik) 



1- 



(27r)2 

v(^')l 



w{k, iJ.;k' , fj,') 



r^ik) 



\Mk)\ 



cos[i?(i;^(fc))-i?(v(fc'))] 
w{k, ii;k' , fj,') 



(23) 



\'^f^<^sm[^{v^{k))-^{v^,{k'))] 



(24) 
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Note that Tjl{k) and T^{k) arc independent of the com- 
mon direction with respect to the angles i9(Vfj^{k)) and 
'd{v^>{k')) in the above integrals are defined, since only 
differences of those angles occur. Therefore, Tj|(fc) and 
{k) are independent of the direction of the electric field 

E. 

Now consider the deviation (7^(fc) of the distribution 
function /^(fc) from equilibrium, 



9f.{k) = Uk) - f. 
Making the ansatz 



(25) 



Mk)\[A^.{k)cos['d{v^{k))] 
+Bi_,{k)sm[div^{k))]^ (26) 

with two parameters vl^(fc), Bn{k) one finds from the 
above equations 



-eErlik) 



1 + 



l(fc) 



B^{k) 



eErlik) 



1 + 



2 ' 



or 



with 



9^.{k)=9l{k)+gi^{k) 



(27) 



(28) 



(29) 



9l{k) = -e\- 



ds 



^E-v^{k), (30) 



df\ r^{k) 



1 + 



{e^ X E^ ■ v^{k) 



(31) 



From this distribution function the electrical current of 
particles in band ji can be obtained as 



3^ = -e 



(fk 



Vn{k)9u.{k) : 



(32) 



and the total electrical current is given by 

From these relations one obtains the following conduc- 
tivity tensor 



xy 



a. 



'xy + CTyy Uyy 

where we have introduced the definitions 

(fk ( df^ 



CTyy — fT, 



± 
XX 

± 



(34) 



(27r)2 



de 



r^{k) 



^k_ 



(35) 



de 



1 



-(fc) 



(36) 



Several remarks are in order: 

(«) For an isotropic dispersion and scattering potentials 



isotropic in real space, only Oxx — o^y are different from 
zero, and the conductivity tensor is proportional to the 
unit matrix. If additionally only one dispersion branch 
is there, the parameter rH becomes 

1 



Tll(fc) 



w{k; k') (l - cos[^{k) - z?(^')]) 



(37) 



This is just the usual expression for the relaxation time 
in the isotropic standard case''® and is independent of the 
wave vector k, T"(fc) =: tq. In fact, the above consider- 
ations can be seen as a generalization of the standard 
isotropic case to general anisotropic dispersions in two 
dimensions. Note that, although the parameters rj and 
have dimension of time, this does not mean that any 
relaxation time approximation has been used to treat the 
case of anisotropic dispersions. In fact Eqs. (29)-(31) 
constitute an exact solution of the Boltzmann equation 
(18). 

(iz) If in addition to aix = f^yy 0, the contributions 
cTxy = cTyx are nonzero, the degeneracy of the conduc- 
tivity eigenvalues is lifted. In the these 
eigenvalues are then given by cr|x ± cr^y with the eigendi- 
rections (1, ±1). 

(Hi) Provided that a^^ = ayy ^ 0, this contribution to 
the conductivity tensor corresponds to the anomalous or 
extraordinary Hall effect. This is an antisymmetric con- 
tribution to the conductivity tensor which does not stem 
from an external magnetic field but entirely from scat- 
tering processes. For such a contribution to be present, 
time reversal symmetry has to be broken. 

For the case of anisotropic dispersions induced by spin- 
orbit coupling as discussed in detail in the next sec- 
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tion, wc will see that there is no anomalous Hall ef- 
fect (since time reversal symmetry is intact), but there 
is a symmetric off-diagonal contribution to the conduc- 
tivity tensor which stems from both aly = crjjx ^ and 



vv 



^0 



IV. CONDUCTIVITY IN THE PRESENCE OF 
SPIN-ORBIT COUPLING 

We now proceed with calculating transport properties 
for Fermi liquid electrons in two dimensions in the pres- 
ence of spin-orbit coupling, using the formalism of the 
previous section. To be specific, we will evaluate the 
transition probabilities in the scattering operator (19) by 
Fermi's golden rule. 



|^|(fc,M|^|fc',M')l' 



(38) 



where V is the operator of a single scatterer and y is 
the density of scatterers. The momentum eigenstates in- 
volved above are normalized as 



{k,ij,\k',n') = A5^p5ij,^ij,. 



(39) 



with A being the area of the system. As a further simpli- 
fication we will consider fixed impurities with (5-function 
shaped scattering potentials (s-wave approximation) , 



V{r) = KS{f) , 



(40) 



where k parametrizes the strength of the potential. The 
square moduli of the matrix elements read 

|(fc,/x|l/|fc',M')l' = Y (l + MM'cos[x(fc) - X{k')]) ■ (41) 

Moreover, we will concentrate on the case of zero tem- 
perature, where the derivative of the Fermi function as 
it arises in the integral expressions for transport parame- 
ters is equal to the negative of a delta function peaked at 
the Fermi energy. Thus the integrations over momentum 
space in Eqs. (35), (36) reduce to integrations over the 
Fermi contour. 

However, even with these simplifications, the integra- 
tions involved are in general non-elementary. In order 
to make analytical progress we concentrate on the case 
of finite Drcsselhaus coupling and small Rashba coupling 
(|a| <C /3, section IV A) and the particular case a = P 
(section IV B). 



A. The case a < /3 

It is straightforward to expand the quantities entering 
the transport parameters and conductivities discussed in 



section HI for a ^ /3 in lowest order in a. However, 
since the calculations are somewhat lengthy, details are 
given in the appendix. The full result for the elements of 
the conductivity tensor (34) up to linear order in a but 
general values for Dresselhaus coefficient (3 and positive 
Fermi energy e/ is stated in Eqs. (A17) and (A18). These 
expressions are still somewhat complicated but simplify 
significantly if one additionally assumes that the "Dres- 
selhaus energy" bd '■= m0^ /It? is small compared to 
as it is usually the case for realistic situations. In other 
words, defining the "Rashba energy" as er := ma^/h^, 
we consider the situation 



£ii < £r> < £/ ; 



(42) 



where the Fermi energy is related to the electron density 
n via (cf. Eq. (12)) 



1 / 2m „ / TO , 



+ O (a^) . (43) 



Then one has 



^xx — ^yy — ^0 



cTxy = cTyx = o-Q (-sign(a)) - 



7 y/eWo 



-1-0 ( £5. [IKi£. 

where uo is the usual Drude conductivity. 



and 



(TO 



e^Tono 



m 



no = 



kj 
2tt 



(44) 



(45) 



(46) 



(47) 



(48) 



are the momentum relaxation time and particle density, 
respectively, in the absence of spin-orbit coupling. The 
eigenvalues of the conductivity tensor are 



— (^xx ~l~ ^xy 



^ — ^xx ^x 



(49) 



with corresponding eigendirections (1, 1) and (1,-1), re- 
spectively. These directions are the symmetry axes of the 
underlying dispersion relations; the same eigendirections 
are found from Eqs. (A17), (A18) where the Dresselhaus 
energy has not been assumed to be small compared to 
the Fermi energy. From Eqs. (44), (45), the conductivity 
anisotropy Act is given by 



Act := 



7 ^srSu 



(50) 
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Wc note that changing the sign of a (by reversing the 
potential gradient across the the quantum well) results in 
a shift by 7r/2 in the wave vector dependence of dispersion 
relations and eigenstates. Such a shift leads to a sign 
change in (cr+ — a~). Therefore, this quantity contains 
only odd powers of a. 

The resistivity tensor p is the inverse of the conductiv- 
ity tensor fulfilling the relation E = pj. From Eqs. (44), 
(45) one finds its components to be 



Pxx — Pyy 

Pxy = Pyx = Posign(a 

+0 



po + C — , — 

7 ^sred 




(51) 



(52) 



with po = 1/co- Thus, a convenient way to experimen- 
tally detect the conductivity anisotropy is a Hall-type 
measurement feeding a current in, say, the a;-direction of 
the quantum well, i.e. = (7xxEx- In the absence of 
Rashba coupling no voltage perpendicular to the current 
is generated, Ey = 0, E^ = Pxxjx- If the Rashba coupling 
is switched on, the off-diagonal elements of the resistiv- 
ity tensor become non-zero, and a finite transverse field 
Ey = Pxyjx = Ex{pxyl Pxx) occurs. We note that this ef- 
fect is similar to the usual Hall effect with the difference 
that the conductivity tensor (and in turn the resistivity 
tensor) is symmetric and not antisymmetric. 



B. The case a = (3 

As discussed already in detail in Ref.^^ and in section 
II, the c:asc; a = fl is special under several aspects. Here 
the transport quantities are readily obtained using the 
form (13) for the dispersion relations. As a result, the 
conductivity tensor is isotropic with 



'vy 



e^rono 



where tq = Ti^ /vn^m as in Eq. (47), and 



\Kt 



no 



27r 



1 2to 



2m/32 



(53) 



(54) 



is the density of electrons. At small deviations from the 
point a = f3 one should expect the conductivity tensor 
to develop again an anisotropy. However, this cannot be 
analyzed in the same way as the case |a| <C /3 since the 
particle velocities and other quantities entering the inte- 
grands in Eqs. (23), (24) do not allow for an expansion 
in |a — /3| around a = (3 for wave vectors with kx = ky. 
At these points the dispersion branches (4) continuously 
merge into two new circles when approaching a = (3, 
cf. the lower right panel of Fig. 1. Therefore, in or- 
der to evaluate the conductivity tensor around a = (3, 



one should use other methods rather than expanding the 
dispersion relations. For our purposes here, we shall be 
content with the statement that the conductivity tensor 
is of course continuous around a = (3, and is isotropic 
exactly at that point. 



C. The influence of trilinear contributions to the 
Dresselhaus term 

The Hamiltonian (2) is derived from the bulk Dressel- 
haus spin-orbit coupling being trilinear in the momentum 

operators^. 



n/bulk 



1 



^ ( cj^'px {pI - pI) + cj^Pv {pI - pI) 
+ ^'Pz{pI-pI)) (55) 



with a coupling parameter 7. In a sufficiently narrow 

quantum well grown along the [001] direction one can 
approximate the operators Pz and p\ by their expecta- 
tion values {pz), (Pz)- This leads to the following two 
contributions to spin orbit coupling resulting from bulk 
inversion asymmetry: The Dresselhaus term (2) linear in 
the momenta with (3 = 7(^2), and the trilinear term 



n 



= l{<T^PxPl-oypypl) 

/(3) 



(56) 



Clearly the typical magnitude of Ti/jj compared to the 
linear term Hd is given by the ratio of the Fermi energy 
£/ of the in-plane motion to the kinetic energy of the 
quantized degree of freedom in the growth direction. For 
typical values of e/ of about lOmeV and not too broad 
quantum wells this ratio is small, and we have therefore 
neglected the Dresselhaus term trilinear in the momen- 
tum components. If desired, it is straightforward to in- 
clude this term in the calculations of transport quantities, 
although the procedure becomes considerably more in- 
volved and will require numerical calculations. However, 
we do not expect, for the following reasons, that including 
the trilinear Dresselhaus term but not the Rashba term 
will lead to anisotropic charge transport: The Hamilto- 



nian 



2m 



'■D 



(57) 



gives the following dispersions for wave vectors k = 
k{cosip, sin (fi): 



21.2 



£±{k,ip) 



h'k 

2m 



±y^/?2fc2 - (1 - cos(4(p)) (^y/37A:4 - y72/ce^ (58) 

The angular variable ip enters only terms of cos(4(^) which 
leads to Fermi contours with fourfold symmetry, dif- 
ferently from the just twofold symmetry in the case of 
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Rashba and linear Drcssclhaus term. In particular, for 
the Hamiltonian (57) the dispersions are symmetric with 
respect to both the axes pairs (1,0), (0,1) and (1,1), 
(1,-1), and these axes pairs are the possible candidates 
for eigendirections of the conductivity tensor. However, 
since directions in the above pairs are equivalent due to 
the existence of the other pair of symmetry axes, we do 
not expect an anisotropy in tranport quantities. More- 
over, as seen above, such anisotropics arise from the in- 
terplay of the Rashba and the Dresselhaus term and are 
tunable by external gates. Concerning gating the Rashba 
coefficient by an electric field across the quantum well, 
one should keep in mind that such an operation might ef- 
fectively also alter the Dresselhaus coefficient /3 = ^{pD 
by changing the expectation value (p^). However, this ad- 
ditional effect cannot change our principle findings con- 
cerning anisotropic transport. 

V. SUMMARY AND DISCUSSION 



understanding of spin-orbit coupling in semiconductors. 

In our calculations we have concentrated on the Dres- 
selhaus contributions being linear in the momentum com- 
ponents, as it is appropriate for not too wide quantum 
wells. The possible influence of the trilinear Dresselhaus 
term is discussed in section IV C. 

The case a = /3 is special imdcr several aspects due 
to the additional conserved quantity that arises at this 
point Here the conductivity tensor is found to be 
isotropic. 

Our approach to anisotropic transport in two dimen- 
sions is based on an exact solution of the Boltzmann 
equation where the drift term is linearized in the in-plane 
electric field driving the current. This formalism can also 
deal with the case of anisotropic single-particle disper- 
sions and should be seen as a generalization of the usual 
isotropic case. We expect this approach to be also use- 
ful in the study of other transport effects such as thermal 
conductivity, magnetothermal effects, and the anomalous 
Hall effect. 



We have presented a theory of anisotropic transport 
in a two-dimensional electron gas. The anisotropy in 
the electrical conductivity is induced by the interplay 
between Rashba and Dresselhaus spin-orbit coupling in 
the scmicondiictor quantum well confining the electron 
gas. The principle axes for anisotropic diffusive charge 
transport are given by the symmetry axes of the single- 
particle dispersion relations, which are anisotropic if both 
Rashba and Dresselhaus spin-orbit interaction is present. 
We have evaluated the conductivity tensor at zero tem- 
perature for scattering on fixed random impurities whose 
potentials are modeled by delta functions. However, be- 
cause of the anisotropic properties of the imderlying dis- 
persions, we do not expect our results to change qual- 
itatively if other impurity potentials are considered. In 
particular, the differential cross section for delta-function 
potentials is isotropic, which makes obivious that our re- 
sult is due to the spin-orbit induced effects and not due 
to special properties of the scattcrcrs. 

To enable analytical progress in the evaluation of 
transport properties, we have concentrated on the case 
of a finite Dresselhaus term and a small Rashba term 
(|a| ^ /3), and on the case where the Rashba and Dres- 
selhaus coefficients are equal (a = (3). For \a\ <C /3 wc 
have found the anisotropic corrections to the conductiv- 
ity tensor due to the presence of the Rashba term. These 
findings demonstrate the principle result that diffusive 
charge transport becomes anisotropic if both Rashba 
and Dresselhaus spin-orbit coupling are present. This 
anisotropy can be tuned by external gates which pro- 
vides the possibility of detecting and investigating spin- 
orbit interaction by measuring spin-unpolarized diffusive 
electrical currents. Apart from possible device applica- 
tions of this effect, the experimental observation of such 
a tunable anisotropy in spin-unpolarized diffusive trans- 
port would certainly significantly confirm and deepen our 
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APPENDIX A: CALCULATION OF TRANSPORT 
PROPERTIES AT \a\ < /3 

In this appendix we present details of the calculation 
of transport properties at \a\ <^ (3 using Fermi's golden 
rule (38) in the case of vanishing temperature. 



1. Dispersion relations, eigenstates, and Fermi 
contours at |a| -C /3 

For the single-particle energies and the phases x(^) 
entering the eigenvectors (5) one finds the following ex- 
pansions: 



X(k) = arg 



-k --k 



a 

I 1 ky ~^ ~p ^ 



7 



= arg ( {-kro + iky) ( 1 - (^x + ikyf 



arg {-k^ + iky) 



exp [ {kx + iky)^] + O 



argi-kx+tky)- - — p 




(A2) 



Here (7^ and ct^ are again usual Pauli matrices acting 
on the two-component vectors k. Note that x(fc) and 
therefore the eigenstates remain unchanged in if k points 
along the directions (1,1) or (1,-1). This can also be 
seen directly from Eq. (6). The expansion of the particle 
velocities reads 



(A3) 



,^ ,7^, , hk B a(k'^a^k) / o\ x 

In order to evaluate the transport parameters tJ. and 
according to Eqs. (23), (24) one needs angles of the type 
t?(v^(fc)). As already remarked in section III, rj. and 
are independent of the direction with respect to which 
these angles are defined. It is convenient to choose this 
direction along the a;-axis. Then one finds analogously 
to Eq. (A2) 




'&{v±{k)) = arg \^(v±{k)^^ + i(v±{k))^ 

2a {P^d^k) 



= 'd{k)± 



O (a2) (A5) 



For the Fermi contour as given in Eq. (10) for positive 
£/ > one has the expansion 



/ 277T, / Tfi \ Tfl 



ma 



T 1 sin {2ip) + O (a^ 



(A6) 



Note that, at a given electron density n, the Fermi energy 

£^ is according to Eq. (43) unchanged in first order in a. 
Moreover, when inserting the Fermi momentum in zeroth 
order in a, 



('i).= 



'2m . / m \ ^„ m _ , , 

Tr^/+ 72 /5'Tt2/3, (A7) 



in (A4) one obtains 



{\v±{kL)\)^ 




(A8) 



i.e. the Fermi velocity is in zeroth order in a independent 
of the band index /z G {+, — }. Using Eq. (A8), the ex- 
pansion (A6) for the Fermi momentum can be rewritten 
up to linear order in a as 



ki {(p;ef) 



{\v±ikL)\)^ 



f3 



It 



. % . sin(2^) I + O {a') . (A9) 
(\v±{kL)\)^ J 



2. Transport quantities 

Using the expansions given in the previous section it is 
a little tedious but straightforward to obtain expressions 
for the transport quantities discussed in section III in 
up to linear order in the Rashba coefficient a. For the 
parameter 1/Tji(k), & {+,—}, one finds in zeroth order 
in a 



1 



rj^^^ 11 + ^- 



/3 







(AlO) 



and the first order is given by 
/ 1 




a0 



(k^a^'k) 



n'2 



{h{\v^ik)\ 



fc2 



(All) 



For 1/T^{k), the zeroth order in a vanishes while the 
first order reads 



1 



I m 



— = VK — 7 — ^ 



2^- 



h 



{\Mk)\)^ 



jk'a'k) 

fc2 



(A12) 



We now turn to the parameters entering the conductivity 
tensor. For the diagonal elements of tr" we find in zeroth 
order 




(A13) 



8 



while the contributions in first order in a vanish. The 
off-diagonal elements of a" are zero at a = 0, and the 
first order reads 

\ ^yj^ V y^Ji h PK^m n n 

-11^ - 21 f §.V 
4 m 16 \h j 




H£/ _L 3 f i\ 
m ~'~ 4 \nj 



(A14) 



Finally, the off-diagonal elements of a-^ vanish up to lin- 
ear order in a, 



+ O (a^) 



(A15) 



The diagonal elements are also zero at vanishing a, and 
the contribution in first order in a is 



h m a (3 j 2e f 



vK?m h h h \ m 




m ^ 16 \h J 



(A16) 



From this one finds the elements of the conductivity ten- 
sor as 



'yy 



e2 



m I 2e 



h vK^m h \ m 



a) h\ m 




(A17) 



h'^ ma 
h vK'^m % h 

2£/ 




(A18) 



Therefore, the eigenvalues of the conductivity tensor are 



xy ) 



'xy 



(A19) 



with corresponding eigendirections (1,1) and (1,-1), re- 
spectively. These directions are the symmetry axes of the 
underlying dispersion relations. 
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FIG. 1. Fermi contours for various values of the Rashba 
coefficient a at a Fermi energy of lOmeV, a Dresselhaus co- 
efficient of lOmcVnm, and a band mass of 0.067 in units of 
the bare electron mass mo. The upper left panel shows the 
isotropic Fermi contours at a = 0. In the upper right and 
lower left panel data at intermediate values of a are plotted, 
while the lower right panel shows data for a = j3. In this case 
the Fermi contours arc two circles having the same radius and 
being displaced from the origin. 
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